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Ground-state fidelity (GSF) and quantum renormalization group theory (QRG) have proven useful tools in the 
study of quantum critical systems. Here we lay out a general, unified formalism of GSF and QRG; specifically, 
we propose a method to calculate GSF through QRG, obviating the need for calculating or approximating ground 
states. This method thus enhances characterization of quantum criticality as well as scaling analysis of relevant 
properties with system size. We illustrate the formalism in the one-dimensional Ising model in a transverse field 
and the anisotropic spin-1 /2 Heisenberg model. 
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Introduction. — Quantum phase transitions (QPTs) take 
place in manybody systems when their ground state or few 
low-lying states of quantum manybody systems — describing 
the systems at zero or almost zero temperature — experience a 
considerable change with variation of Hamiltonian parameters 
Jjl. The standard symmetry -breaking mechanism most often 
fails to capture QPTs. In fact, it is not always clear how to de- 
fine a suitable local "order parameter" signifying a symmetry 
breaking at a critical point; e.g., to systems exhibiting "topo- 
logical order," no local order parameter can be attributed |2|]. 
Moreover, discontinuities or singularities of the ground-state 
energy cannot always predict QPTs 0]. 

Such inherent difficulties with identifying QPTs have ne- 
cessitated tools that could better capture nature of quantum 
correlations. Along these lines, "entanglement" has proved 
a useful signature for some QPTs J3]. More interestingly, 
though, the elementary concept of the "ground-state fidelity" 
(GSF) has recently been shown to provide another remarkably 
useful means in signaling QPTs 15J. This may be somehow 
natural as the ground state encodes all relevant information 
about a quantum system at zero temperature, hence a phase 
transition is expected to be identified by, e.g., a considerable 
difference between the ground states right before and right af- 
ter a quantum critical point. This enables GSF as a fairly gen- 
eral order parameter for quantum critical systems (irrespective 
of their internal symmetries) l5j,|6J], endowing as well a rich 
intrinsic geometric feature 121 [8[] . 

Alternatively, "quantum renormalization group" (QRG), a 
variant of RG at zero temperature @], puts forward a tractable 
recipe for studying critical behavior of a variety ofquan- 
tum manybody systems, especially in one dimension I loT 121. 
QRG essentially hinges on a coarse-graining procedure un- 
der transformation of Hamiltonian parameters, to weed out 
irrelevant short-distance information while retaining original 
large-scale picture after rescaling length. This formalism has 
recently been employed successfully to find critical properties 
of a variety of quantum manybody systems |l3 - 151. 

Despite the utility of GSF, in practice its applicability is 
largely restricted to systems for which one can somehow com- 
pute ground state or an approximation thereof — which is a 
demanding task. To overcome this issue with computation 
of GSF, approaches based on, e.g., tensor network 11611 and 



Monte Carlo algorithms [ 17] have recently been employed. 

Here we observe that QRG can also offer a powerful al- 
ternative approach to computing GSF. Specifically, we aim to 
combine the GSF and QRG formalisms into a unified picture 
for identification of QPTs, obviating the need for the knowl- 
edge of ground state. Our formalism is fairly general and in 
principle can be applied to a broad class of manybody sys- 
tems which are amenable to QRG formalism. To illustrate the 
framework, we elaborate it within two examples: (i) the Ising 
model in transverse field (ITF) and (ii) the anisotropic (XXZ) 
Heisenberg model. For specificity, in the following we adopt 
the Kadanoff RG recipe |9|], although the formalism is appli- 
cable to other RG schemes as well. 

Formalism. — Consider a quantum system of N spins (each 
with the Hilbert space W s of dimension s), defined on a 
Hilbert space "Hf = Hf N , with the Hamiltonian H(x), 
where x (for simplicity taken to be a single parameter) is a 
coupling constant. In the renormalization procedure, the orig- 
inal model Hamiltonian H is replaced with an effective or 
renormalized Hamiltonian H at the cost of renormalizing cou- 
pling constants JlolfTl . 131. As a result, the original Hilbert 
space is also mapped into a renormalized Hilbert space 
H encompassing only the effective degrees of freedom. Inte- 
grating out less important degrees of freedom gives rise to a 
flow in the coupling constant space. We can define an em- 
bedding operator T(x) : W — s- H to represent this step: 
T(x)|$«(x)) 



|$o(x)), where |$ ) and |$o^) are the 



ground states of H and H = T T HT, respectively. T is usually 
constructed as follows: divide the system lattice into blocks of 
a given size, say, m; considering the original Hamiltonian, at- 
tribute a Hamiltonian hf to each block /; diagonalize hf to 
find eigenvectors {\4>i)i}f =1 corresponding to first s eigenval- 
ues to formT/ = Yh=i \4>i) I (4>i\, where {|0i)j}f =1 represent 
new block degrees of freedom constituting Hg (and overall 
H^ m ); and, finally define the global embedding operator as 



T 



We now recall the definition of "fidelity" /, for a system of 
size N < oo, associated to the ground states |$q( x ±)) as 

f = f(x,5;N) = ($ (x-)|$o(x+)) (1) 
= <4 1) (x_)|Tt(x_)T(x + )|< P ( 1) (x + )), (2) 
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where x± = x ± S, and 8 represents a small varia- 
tion of x — dropping its customary absolute value for now. 
The group property of the renormalization procedure en- 
sures that ^(x)) = |$ (x (1) )), where xW = xW(x) 
is the renormalized coupling. Hence, the fidelity can be 
written in terms of the renormalized coupling as / = 
($o(x a) )|Tt(x_)T(x + )|$ (x ( f 1) )). In some cases, the 
right-hand side may be written as a function of /W, leading 
via RG iterations to a recurrence relation of the generic form 
= (£ > 0), where R is a model-dependent 

function, and f^> is the GSF after £ RG iterations. Solving 
this equation (analytically if the model is amenable to some 
exact methods) and utilizing a priori knowledge of the associ- 
ated QRG fixed points x c can provide useful information such 
as behavior of the GSF in/around a quantum critical point or 
how it scales with the system size. 

Equation (O indicates that a significant simplicity ensues 
for the cases in which T t (x^)T(x + ) = w (0) (x, (5)1 with 
some u/°) (hence R is linear) 11811 : specifically, here / = 
w (o)y(i) The RG iteration yields 



/(x,<5;iV) = /W(x, ( 5)n^ ' ) (x, ( 5), (3) 



3=0 



reducing the computation of the GSF to and 



in which N = 



U+i 



An immediate consequence is 
that if |u>W(x,<J)| < 1 W, then lim^^ /(x, S) = 
[note lim^oo = limjv-)-oo]> implying an "orthogonality 
catastrophe" — hence QPT — for the corresponding x. If for 
a continuum of xs such a behavior persists, we have a critical 
line (as in the XXZ model discussed later). 

It is often the case that rather than the GSF, the GSF 
"susceptibility," defined through Taylor expanding / up to 
0(S 2 ), f(x,6;N) ps 1 - (£72)x(x; N), suffices to cap- 
ture quantum criticality [0, U^. Expanding the embedding 
operator and using the identity Tt(x)T(x) = 1 yield / ps 



($ (xL 1) )|5(x,5;iV)|$ (x! f 1) )),where5 



-2(5Ttd x T- 



25 2 d x T^d x T. Thus the above RG procedure for / applies to 
X as well provided that we can treat S appropriately. 

The ITF model. — This model on a periodic chain of N sites 
is defined with the Hamiltonian 



H(J,g) 



N 



(4) 



where J defines an energy scale, g is the parameter that con- 
trols QPT, and af is the Pauli matrix for site i. To apply QRG, 
the chain is divided into blocks of m = 2 sites described by 
h b = Y^jllhf, where hf = - JK^afj + gofj) O. 
The block-block interaction Hamiltonian is also represented 

by H BB = -JEfi?( cr I,/ cr i,/+i+9 cr 2j)- h f canbediago- 
nalized exactly, whence Tj is constructed from the two lowest 
eigenstates as Tj = |</>i)/(fr | + \<t>2)i(fy \, in which |^ x )j = 
A(g)\ ft) + B{g)\ U), \^)i = A(g)\ ft) + B(g)\ |t> are 
the two degenerate ground states of hf. Here {| f), | \)} 



are the eigenvectors of a x , {\ ^}i,\ JJ-)/} represent the states 
of block I, and A(g) = s/y/l + s 2 ', B(g) = l/Vl + .s 2 ', 
with s = g + yl + g 2 . The global embedding operator 
T = ®^T T leads to H = T^(H B + H BB )T, which is 



akin to the original one [Eq. (0]i] modulo replacing the cou- 
pling constants with 



JW = J- 



2(g+y/T+F) 
1 + (.9 + Vl+F) 2 



g (1) = g 2 - 



(5) 



The renormalized couplings after I RG iterations are obtained 
simply from Eq. @ upon substituting (J, g, J^ 1 ', g^) — > 
(j(t-i) i g(t-i) ) We note that for this model the 

RG fixed points are 6 {0, 1, oo}, from which g c = 1 is unsta- 
ble whereas and oo are stable under the RG flow: 



» « * » 







1 



oo 



A straightforward algebra shows that here 

fW(g,6) = A(g W )A(g^) + B{g^)B(gf) and 

^(g,S) = [A{g^)A{gf) + B{g^)B{gf)] N/23+ \ 
from whence through Eq. (0 one can find an analytical 
expression for fnp(g,S; N). Figure Q] shows fup(g,S; N) 
for various values of (5, N). A drop is seen at g = 1, which 
verifies it as a quantum critical point. In agreement with 
Ref. i20ll . two regimes NS < 1 and NS > 1, corresponding 
respectively to the "small-size limit" and the "large-size 
limit," can be discerned. In the small-size limit, the GSF 
drops at/around g c = 1; whereas, in the large-size limit, the 
GSF drops to zero for a domain of gs around g c = 1 whose 
size is ps 0(5) (Fig. Q] [right panel]). 

In the thermodynamic limit, an interesting limit of the GSF 
can also be obtained. Equation ([5]) dictates three distinct re- 
gions for g± before applying the 6 — > limit: (i) g_ < 1 and 
g+ < 1, (ii) g- < 1 and g + = 1, or g_ = 1 and g + > 1, (iii) 
g_ > 1 and g + > 1. Thus considering the RG flow diagram 
of the ITF model leads to 



lim lim f nP (g,8;N) 

5— s-0 AT->oo 



1 



1 + 
V4 + 2V2 1 ' 



1 <W6) 



(where 8 g \ is the Kronecker delta function) which seems to be 
consistent with the small-size limit (Fig. Q] [left panel]). Note 
that this sharp drop of the GSF is a signature that g c = 1 is a 
quantum critical point, consistent with what QRG suggests. 

Moreover, the drop in the GSF accordingly signals a non- 
analyticity in Xitf, implying that here the GSF susceptibility 
is also a reliable to identify the criticality of the model. In 
fact, expanding f(g,6;N) up to 0(S 2 ) (as explained in the 
formalism) yields the following expression: 



x [l) {g) 

2 e+i 



(d e A(gW; 



diB(g 



(O 



n 



(7) 



with d ( = dgw, thence Xrrp(g) = X (i) {g) + 
NEU D(gW)/2?-\ where D(g) - s 2 /[(l+ 3 2 )(l+ s 2 ) 2 ]. 
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FIG. 1. GSF of the ITF model vs. the field strength g. Left: Small-size limit NS < 1. Right: Large-size limit N5 > 1. 



{t-l) 



fir 



1, = N 2 /32, hence 



Note that = 

Xitf(<?c) ~ N 2 |fj. Alternatively, for second-order QPTs, 
it has been known that x(ffc) ~ N 2 / dv ', where d is the di- 
mensionality and v > is the critical exponent capturin g th e 
divergence of the correlation length £ ~ \g — g c \ ~ v j8l l2lll22ll . 
Thus we obtain i^tf = 1 . 

Further scaling analysis can be made for /itf(<7, 8; N) with 
fixed N or 8. Figures|2]and[3]show ln(— In /) close to g c = 1 
vs., respectively, In 8 for fixed N = 32768 and In N for fixed 
8 = 0.001, comparing two ground states with gx = 1 — kS 
and g 2 = 1 - (fc + 1)<5, where fc e {0, 1, 2, 3, 4}. In Fig. El 
the k — case, labeled by (1, 1 — 8), comparing the ground 
state at the quantum critical point and a state very close to it in 
the ferromagnetic phase, shows a behavior akin to the k ^ 
cases for small 8s in the N8 < 1 regime; however, it shows a 
distinct behavior for N8 > 1, signaling that one of the states 
is exactly at the quantum critical point. The k ^ cases, 
comparing in fact two states in the ferromagnetic phase close 
to the quantum critical point, exhibit the In / ~ — 8 2 scaling 
for N8 < 1, connected with a crossover domain of N8 as 
0(1) to the In/ ~ —8 scaling for N8 > 1, in agreement 
with Ref. Q. In contrast, Fig.[3]exhibits the In / N 2m 



scaling for N8 < 1, connected with a crossover domain of 
N8 as 0(1) to the different lnf ~ —N scaling for N8 > 1, 
also in agreement with Ref. 12011 . The distinct behavior of the 
(1,1 — 8) case hints that one of the states is at the quantum 
critical point. The scaling behavior of the GSF around g c = 1 
in the paramagnetic phase (g > 1) gives results similar to 
Figs. [2] and [3] which we did not present here. 

The XXZ model. — The anisotropic spin-1/2 Heisenberg 
model (XXZ) on an open chain is defined with 12311 



N 



H(J,A) = jJ2h 



(8) 



where J > is the exchange-energy coupling, A = (q + 
q~ 1 )/2 is the axial anisotropy given in terms of a pure phase q, 
and/i M+ i = afaf +1 +af<Tf +1 +a + {q)a*tTf +1 -a-(q){a* - 
°f+i)> w i tn a ±(l) — (l i 9~ 1 )/2- Equation (O is different 
from the ordinary XXZ Hamiltonian in the boundary term cx 
°i ~ °fr> which is unimportant in the thermodynamic limit. 
It is known that this model is critical (gapless) for |A| < 1 
(critical line), exhibiting no long-range order |l]]. 

The RG procedure here is implemented based on the quan- 




ln(5) ln(N) 



FIG. 2. Scaling behavior of the GSF for the ITF chain with fixed size 
N = 32768: ln(- In /) vs. In S close to g c = 1. 



FIG. 3. Scaling of the GSF for the ITF chain with fixed 5 = 0.001: 
ln(— In /) vs. In N close to g c — 1. 
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FIG. 4. Average GSF susceptibility vs. 1/N for the XXZ chain. 
In the inset the horizontal axis is in the normal scale (unlike the log 
scale of the main plot). The behavior in the inset is reminiscent of 
Fig. 1 of Ref. I2II1 obtained through an exact diagonalization. 



turn group property of the Pauli matrices 12311 . The Hamil- 
tonian ([8]) is decomposed to 3-site blocks (m = 3), where 
block / is comprised of sites interacting as hf = /i^.^+i + 
hi I+ i^ I+ 2, and the rest of the Hamiltonian constitutes the 
block-block interaction. The ground state of the block Hamil- 
tonian is doubly degenerate, represented in the a z -basis as 



\<h) = -c+\- 

\fo) 



-c 4 



C-) 



+)-c_\- 

C-\- 



-+>> 
-+>, 



The embedding opera- 



where C± = yjq^/iq + q' 1 + 4) ', 
tor for block / is then given by Tj = 
where {| ft)/, | JJ-)/} denote states of block /. Therefore, the 
renormalized Hamiltonian is obtained similar to Eq. (0 with 
the following renormalized coupling constants: 



q 



AW = A , J (1) = J 



A straightforward calculation shows that here 



'(A, 5) 



2V(2 + A_)(2 + A + )' 



•(9) 



(10) 



where 



u(A, S) = V2 2Jl + A_A+ + J {I - A 2 _)(l - A\ ) 



1 + A_A + - X /(1-A 2 _)(1-A 2 



.(11) 



Now through the RG formalism (especially noting that A^ = 
A), the following analytical expression is obtained for the 
GSF given any (A,5;N): 



(12) 



Figure HJrepresents this scaling for some values of A and 5. 

Interestingly, this simple and elegant relation also enables 
detection of the associated criticality in the XXZ model (wit- 
nessed in Refs. 12 ll 12411 ). resolving conclusively the doubt 



that the GSF might be insufficient 12511 . Evidently we have 
cjW = and \ojW | < 1 for all S ^ and |A| < 1; hence 

lim Urn fxxz(A,5;N) = 0, V|A| < 1. (13) 

(5— s-0 N->oo 

That is, the whole |A| < 1 line is critical HI]. This is a re- 
markable result in that to characterize the criticality of the 
XXZ model we did not need to know the ground state or an 
approximation of that j2lll24Tl . We remark that the GSF sus- 
ceptibility xxxz = 6(7V - 1)/[(1 - A 2 )(2A + 4) 2 ], obtained 
through our RG approach, in contrast to the GSF only cap- 
tures the criticality at the symmetric point | A| = 1 (including 
the Kosterlitz — Thouless point A = 1 and the ferromagnetic 
critical point A = — 1). 

Summary. — We have developed a viable, general quantum 
renormalization group formalism to calculate ground-state fi- 
delity in quantum critical systems. Our formalism combines 
two powerful methods in a unified framework, enhancing 
characterization of criticality in quantum manybody systems. 
Specifically, our formalism is structured on coarse-graining a 
quantum system (e.g., by partitioning it into blocks) and then 
rescaling system length in order to eliminate short-scale or ir- 
relevant interactions from Hamiltonian. In various cases this 
enables a renormalization-based recurrence relation for the fi- 
delity, without the need to know system's ground state, an 
approximation thereof, or an order parameter. With this ad- 
vantage, one can utilize the quantum renormalization group 
toolkit to boost or even simplify calculation of critical proper- 
ties in systems where renormalization works sufficiently well. 

We have illustrated our formalism through two examples, 
the Ising model in transverse field and the anisotropic Heisen- 
berg chain. In both models, our approach produced analytical 
expressions for the ground-state fidelity, resulting to the ex- 
pected criticality (especially in a simpler and more conclusive 
way than already had been suggested for the second model). 
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